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ABSTRACT 

The flare/spotted spectroscopic binary star CC Eri was observed with the Position 
Sensitive Proportional Counter (PSPC) on the X-ray satellite ROSAT on 1990 July 
9-11 and 1992 January 26-27. During the observations, the source was variable on time 
scales from a few minutes to several hours, with the X-ray (0.2-2 keV) luminosity in 
the range ~ 2.5 — 6.8 x 10 29 erg s _1 . An X-ray flare-like event, which has a one hour 
characteristic rise time and a two hour decay time, was observed from CC Eri on 
1990 July 10 16:14-21:34 (UT). The X-ray spectrum of the source can be described by 
current thermal plasma codes with two temperature components or with a continuous 
temperature distribution. The spectral results show that plasma at T e ~ 10 7 K exists 
in the corona of CC Eri. 

The variations in the observed source flux and spectra can be reproduced by a flare, 
adopting a magnetic reconnection model. Comparisons with an unheated model, late 
in the flare, suggest that the area and volume of the flare are substantially larger than 
in a solar two ribbon flare, while the electron pressure is similar. The emission measure 
and temperature of the non-flaring emission, interpreted as the average corona, lead to 
an electron pressure similar to that in a well-developed solar active region. Rotational 
modulation of a spot related active region requires an unphysically large X-ray flux in 
a concentrated area. 
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Ol INTRODUCTION 

Kfl CC Eri (HD 16157) is a spectroscopic binary with period 
^ of 1.56 days. It consists of a K7Ve (although some authors 
adopt MOVe) primary and a ~dM4 secondary with mass 
ratio PS 2 (Evans 1959, 1971; Strassmeier et al. 1993). The 
primary co-rotates with the orbital motion owing to the tidal 
lock and is one of the fastest rotating solar neighbourhood 
K dwarfs known. The stellar parameters adopted for the pri- 
mary are given in Table 1. (See Section 4.1 for the definition 
of Ro, the Rossby number). 

CC Eri is classed as a very active BY Draconis variable 
flare star (Evans 1959; Busko & Torres 1976). The system 
shows H Balmer and Ca II H and K lines in emission (Evans 
1959). It exhibits sinusoidal optical light variations of vary- 
ing amplitude and phase with a periodicity essentially iden- 
tical to that of the orbital motion (Evans 1959). The varia- 
tions have been attributed to the rotational modulation of 
dark star spots (Krzeminski 1969; Bopp & Evans 1973). The 
optical emission line fluxes vary roughly in anti-phase with 
its photometric variation (Busko, Quast & Torres 1977). The 
source is also visible in the far-infrared wave-band with a lu- 
minosity at 12 /tin of ~ 0.22 x 10 31 erg s" 1 (Tsikoudi 1988). 
Radio observations indicate that CC Eri is variable at ra- 



dio wavelengths and is in a strong flaring state for certain 
periods of time (Caillault, Drake & Florkowski 1988; Giidel 
1992). 

CC Eri has also been observed at ultra-violet wave- 
lengths with IUE and some evidence of flaring was noted 
(Byrne et al. 1992). A multi-wavelength campaign in op- 
tical and ultra-violet has shown that CC Eri is the most 
luminous solar neighbourhood K star yet measured in the 
C IV lines (Byrne et al. 1992). The observed variation of 
the Mg II lines with phase confirms earlier conclusions that 
the chromospheric emission is anti-correlated in a general 
sense with the optical continuum, suggesting the presence 
of active region emission associated with star spots. 

Previous X-ray observations of CC Eri have been 
made with HEAO-1, Einstein <m& EXOSAT (Tsikoudi 1982; 
Schmitt et al. 1987, 1990; Pallavicini et al. 1988; Pallavicini, 
Tagliaferri & Stella 1990). These observations yield a source 
luminosity between 1.6 and 4.2 x 10 29 erg s -1 and temper- 
atures of 4.6 to 15xl0 6 K. Since none of these observations 
lasted very long, little information on the temporal and spec- 
tral variation of the source was available. 

Taking advantage of the increased sensitivity and the 
low background noise of the ROSAT PSPC detector, com- 
pared with Einstein and EXOSAT, we have obtained high 
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Table 1. Adopted parameters of CC Eri. 

Spectral Type d R*/Rq log P R T e s 

(pc) (CGS) (days) (ergs" 1 ) (K) 

K7Ve+dM4 12 0.60 4.55 1.56 6 x 10" 2 2.80 x 10 32 3875 



quality, low resolution spectra of CC Eri. In this paper, we 
present the first X-ray observations that reveal details about 
the spectrum of the source and its time variations. We de- 
scribe the PSPC observations of CC Eri and the data anal- 
ysis methods in Section 2 and present the observational re- 
sults in Section 3. In Section 4 we show that the variations 
in the observed flux can be interpreted in terms of a flare 
using magnetic reconnection theory. In Section 5 we use the 
emission measures and temperatures to make an alternative 
model of the flare decay phase, assuming that no heating 
occurs, and also model the quiescent parameters. Finally 
in Section 6 we summarize the results, and suggest further 
observations. 



2 OBSERVATIONS AND DATA ANALYSIS 

The ROSAT PSPC observations of CC Eri were made on 
1990 July 9-11 during the ROSAT calibration and verifica- 
tion phase, and on 1992 January 26-27 on our behalf during 
the Guest Observer Program AO-1 phase. The Boron filter 
was used in the observation made on 1992 January 26-27 
to allow an increase of spectral resolution at lower energies. 
The total exposure time was 17718 seconds, divided into 22 
groups for the study of spectral variability discussed below. 
Table 2 lists details of the observations and the time inter- 
vals over which the X-ray spectrum was accumulated. The 
PSPC has an energy range of 0.1-2.4 keV with a spectral 
resolution of AE / E Rj 0.42 at 1 keV, more than twice that 
of the Einstein Image Proportional Counter (IPC). A com- 
prehensive description of the ROSAT satellite, the X-ray 
telescope and the PSPC detector can be found in Trumper 
(1983) and Pfeffermann et al. (1987). 

The observational data were reduced using the Star- 
link ASTERIX X-ray data reduction package (Version 1.6b). 
The X-ray counts of the source were extracted from a cir- 
cle, centred on the source with a radius chosen to optimize 
the signal-to-noise ratio (6 arcmin in the case of CC Eri). 
Background counts, estimated from an annulus region cen- 
tred on the source cell with a larger radius (15 arcmin), 
but without obvious contaminating sources, were subtracted 
from the source counts. Standard corrections for shadowing 
by the wire mesh, mirror vignetting, photon scattering and 
also dead time were applied to the background subtracted 
source counts. 

The X-ray pulse height spectra were analyzed using 
the XANADU-XSPEC package (version 8). After exclud- 
ing "bad" raw pulse height analysis (PHA) channels cor- 
responding to energies below 0.1 keV and above 2.4 keV, 
the 256 channel data were binned up according to the com- 
pressed channel assignment of the ROSAT Standard Anal- 
ysis Software System (SASS). The rebinned spectrum con- 
tains 34 SASS channels and each new bin has a constant 
oversampling factor (2.45) of the spectral resolution. The 
SASS channels were selected to prevent dramatically dif- 
ferent oversampling factors at various energies skewing the 



Table 2. CC Eri - ROSAT observation journal. 
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spectral fit. The latest recent version (released in 1993 Jan- 
uary) of the ROSAT PSPC response matrix was used in our 
spectral analyses. A systematic error of one percent of the 
measured count rate, in each of the 34 SASS channels, was 
quadratically added to the photon noise in those channels 
to allow for uncertainties in the calibration of the PSPC 
response matrix. We further restrict our data analysis in 
the SASS channel 5-30 (~ 0.17 — 2 keV) to avoid possible 
soft X-ray contamination (solar-scattered X-rays, auroral X- 
rays, etc) below 0.17 keV and the large uncertainty in the 
calibration (systematic error > 10%) above 2 keV. 

Various X-ray plasma emission codes are available for 
modelling the X-ray spectra of late-type stars. We have 
adopted the most recently available versions of two of them, 
one by Landini and Monsignori-Fossi (LMF) and the other 
by Raymond and Smith (RS), to model our data in terms 
of the emission measures and temperatures. Discussions of 
these two models can be found in Landini & Monsignori- 
Fossi (1970, 1985, 1990) and in Raymond, Cox & Smith 
(1976), Raymond & Smith (1977) and Raymond (1988). 
Both models assume an optically thin plasma in statistical 
equilibrium and include emission from important elements 
in addition to continuum radiation from bremsstrahlung, re- 
combination and two-photon processes. Given the low spec- 
tral resolution of our PSPC data we have frozen elemental 
abundances at values given by Anders & Grevesse (1989) 
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and Grevesse, Noels & Sauval (1992) for the solar photo- 
sphere. We have not attempted to model the spectral lines 
in detail. 

When applying the LMF and RS codes to fit our data, 
we assume that the X-ray emission from the star can be 
described by the form 



F(T e , E) 



1 



4tt(2 2 



Nl&VP(T e ,E) erg cm" 2 s" 1 keV _1 (l) 



where d is the source distance, N e is the electron density, and 
P(T e , E) is the plasma emissivity as given by the LMF or RS 
codes at electron temperature T e and energy E. The volume 
emission measure, Em(V) = J N^dV , derived, assumes that 
all photons escape, none being intercepted by the presence of 
the star. This is referred to later as the "apparent" emission 
measure. 

We have used several different combinations of the RS 
and LMF models to fit the PSPC spectra to derive tempera- 
tures and emission measures by minimizing the \ 2 statistic. 
These combinations include one, two, and multiple temper- 
ature models (i.e. the X-ray emission is F(T t , E), n = 
1,2,3...), and a continuous emission measure distribution 
(based on the RS plasma code) in which the volume emis- 
sion measure at a given temperature T e is defined as 
Em(V,T e ) = £m(TOmax(T e /T max ) a , where Em{V) m ^ is 
the volume emission measure at the maximum temperature 
Tmax, and a is the gradient. Using the differential emission 
measure dEm(V, T e )/dT e = NgdV/dT e , the X-ray flux from 
the star in the continuous emission measure distribution 
model can be expressed as 



-F(Tmax, E) 



aEm(V) n 



1 



4tt(2 2 



J- m 



T™' 1 P(T e , E)dT e erg cm" 2 s" 1 keV" 1 (2) 

T min 

where T m ; n is the minimum temperature. In the Sun the 
emission measure distribution as a function of T e has a 
roughly constant gradient down to T e ~ 2 x 10 5 K, so this has 
been adopted for T m ; n . However, the spectral range included 
by the PSPC contains little emission from lines formed be- 
low 4.7 x 10 5 K, and the consequences of using this higher 
T m ; n are discussed in Section 3.2. 

It is important to stress that the one, two, and multiple 
temperature models provide only a parameterization of the 
coronal emission and do not necessarily represent a physi- 
cally consistent description of the physical processes in the 
emitting plasma. Observations of the quiet Sun and solar 
active regions show that a continuous emission measure dis- 
tribution is present (e.g. Malinovsky & Heroux 1973; Dere 
1982). 

We find that the one-temperature model gives an un- 
acceptable fit to all of the PSPC spectra and that including 
more than two temperatures does not improve the statistics 
of the spectral fit. Therefore these two types of models will 
not be discussed any further in this paper except when pre- 
senting the EinsteinYPC results where the IT model gives an 
adequate description to the low spectral resolution IPC data. 
In general the two-temperature (2T) RS and LMF model 
and the continuous emission measure (CEM) RS model give 
equally good fits to the PSPC spectra. 

Interstellar absorption is considered using the model of 



Morrison & McCammon (1983). However, the hydrogen col- 
umn density Nu , when allowed to vary as a free parameter in 
the fit for the 2T and CEM models, is not well constrained 
by the compressed 0.17-2 keV PSPC spectrum. We have 
therefore fitted the 2T and CEM model to the data with 
column density N H fixed at 0.0, 2.6, 10, and 29 xlO 18 cm" 2 . 
All the models give an acceptable fit to the data. We note 
that varying Nu in the range of 0.0-29 xlO 18 cm" 2 has lit- 
tle or no effect on the temperature and emission measure of 
the high temperature component of the 2T model. However 
it causes an average change of ~ 10% in the temperature 
(which decreases as TVh increases) and in the emission mea- 
sure (which increases as TVh increases) of the low temper- 
ature component. These small changes, in most cases, are 
within the one-sigma uncertainties of the temperature and 
emission measure (see Table 3). For the CEM model, both 
the gradient a and the emission measure are correlated with 
the change of the low energy absorption and the temperature 
Tmax does not change with Nu- Since Nu = 2.6 x 10 18 cm" 2 
is estimated from an average interstellar hydrogen density of 
0.07 cm" 3 (Paresce 1984) and the distance of CC Eri, in the 
following sections we will only present the spectral results 
derived with this column density. 



3 RESULTS 

3.1 X-ray Intensity Variability 

CC Eri was variable on time scales from a few minutes 
to several hours during the PSPC observations. The X- 
ray (0.2-2 keV) luminosity of the source was in the range 
~ 2.5 — 6.8 x 10 29 erg s" 1 . We have re-analyzed the X-ray 
observations of CC Eri made with the Einstein IPC and the 
EXOSAT Low Energy (LE) telescope in order to compare 
our ROSAT results with them in a similar energy range. 
The IPC measurement gives a luminosity of ~ 2.6 x 10 29 
erg s" 1 in the 0.2-2 keV band and the EXOSAT LE gives 
~ 6.5 x 10 29 erg s" 1 in the same energy band, both of which 
lie within the range of values from the ROSAT observations. 

Fig. 1 shows the X-ray light curve of CC Eri obtained 
on 1990 July 10 16:14-21:34 (UT). The data are in seg- 
ments separated either by earth occultation of the source or 
portions of orbits with high background. A flare-like event 
was observed in this period. The source became gradually 
brighter and reached a peak flux of 6.3 ±0.1 counts s" 1 (cor- 
responding to a luminosity ~ 6.5 x 10 29 erg s" 1 ) at around 
17:53 (UT). The source flux slowly decreased again and after 
about 3.7 hours the luminosity became ~ 3.3 x 10 29 erg s" 1 . 
There is a bump in the light curve between 19:39-19:57. If 
the source behaviour did not change dramatically during the 
gaps, the observed flare-like event has a one hour (or less) 
characteristic rise time and a two hour decay time. The ratio 
of the peak to minimum fluxes is about 2. 



3.2 X-ray Spectrum Variability 

Figs. 2(a) and 2(b) show respectively the two X-ray spec- 
tra of CC Eri, obtained on 1990 July 10 17:47:29-18:00:49 
(spectral group 7) when the source was at the peak lumi- 
nosity, and on 1992 January 26 21:08:31-21:38:31 (spectral 
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Figure 1. X-ray light curve of CC Eri on 1990 July 10. Each data point represents the (0.2-2 keV) X-ray flux (counts s — 1 ) in a 100 
second time bin. The one sigma error is given for each data point. 



group 21). The boron filter was used in obtaining the spec- 
trum shown in Fig. 2(b). The spectra are shown in terms of 
the detector counts (crosses). The fits with the 2T RS model 
are shown as histograms. The residuals of the model are also 
given in Fig. 2 in terms of the significance which is defined 
as (data— model)/one-sigma error. It is obvious from Fig. 2 
that the 2T RS model gives a good fit to the spectra. The 
fits using the CEM RS code are not shown, but are virtually 
indistinguishable in the fits and residuals. The parameters 
of all 22 PSPC spectra are listed in Table 3 for the 2T RS 
and LMF models and in Table 4 for the CEM RS model. Ti , 
Em(V)i, T2 and Em(V)2 in Table 3 are the temperatures 
and volume emission measures of the first and second com- 
ponents. The spectral results from fitting the Einstein IPC 
spectrum with the IT RS & LMF models and the CEM RS 
model are also listed in Table 3-4 for the purpose of com- 
parison. The errors quoted in both the tables and figures 
discussed below correspond to a 68% confidence level. 

Since the available codes give different results we have 
examined the relation between the model parameters pro- 
duced. As can be seen from the numbers given in Table 3 
there are systematic differences between the spectral results 
derived with the RS and LMF plasma codes although they 
both give adequate fits to the spectra and give a similar 
relationship between the source spectral parameters and its 
flux. The emission measures of the high temperature compo- 
nent are correlated with the source flux. Because the LMF 
code leads to larger high temperature emission measures, 
the smaller, low temperature component is less systematic 
in its behaviour. The temperatures derived from both the 
RS and LMF models cover a smaller range of values than 
do the emission measures. The variation of T2 from the RS 
code (but not from the LMF code) increases slightly with 
the X-ray flux. The temperature T\ does not vary systemat- 



ically with the flux and has an average level of ~ 2.3 x 10 6 
K for the RS model and with ~ 1.8 x 10 6 K for the LMF 
model. The IPC temperatures and emission measures given 
in Table 3 are within the range of the PSPC measurements. 

There is a clear correlation between the ratio of the high 
temperature emission measures and the temperature T2 , as 
shown in Figure 3. Em(V) 2 (LMF) / Em(V) 2 (RS) decreases 
systematically as either T2(LMF) or T2(RS) increases. The 
situation regarding the low temperature component, how- 
ever, is more complex. In most cases for the same PSPC 
spectrum when the LMF model gives a lower temperature it 
will at the same time yield a lower emission measure than the 
RS model. These systematic differences demonstrate the sys- 
tematic, temperature dependent uncertainties in the plasma 
emissivity codes currently available, e.g. in the ion popula- 
tions adopted. 

The numbers given in Table 4 show that the plasma 
temperature T max of the CEM RS model is in the range of 
11-32X10 6 K and is roughly correlated with the source flux. 
The emission measure £m(F) miI shows a small increase 
with the flux (Note that the two highest values of £m(F) miI 
come from spectra obtained using the Boron filter). The 
gradient a, which does not correlate with the flux, is spread 
around a weighted average value of 1.15 ± 0.02. We have 
fitted the Einstein IPC spectrum using the CEM RS model 
with Nh fixed at 2.6 x 10 18 cm -2 and a at 1.15. As given 
in Table 4 these measurements of T max and £m(V) miI are 
consistent with the PSPC results. The spectral results from 
the CEM RS model depend on the value of the minimum 
temperature T m ; n . Changing T m ; n from 2 x 10 5 K to 4.7 x 
10 5 K does not change T max and Nu, but would increase 
_Em(y) max by ~ 40% and a by ~ 10%. Given the present 
uncertainties in the codes we do not at this stage attempt 
to interpret the value of a. 
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Figure 2. X-ray spectra of CC Eri fitted with the 2T RS model, (a) Data obtained on 1990 July 10 17:47:29-18:00:49; (b) Observation 
on 1992 January 26 made with the Boron filter. The top panels show the spectra in terms of the detector counts (crosses, in units of 
counts s — 1 keV - 1 ), and the 2T model (solid histogram). The residuals of the model fitting are given at the bottom of (a) & (b). 




Figure 3. Correlation between Em(V) 2 (LMF) / Em(V) 2 (RS) and X 2 (LMF) (filled triangles) or T 2 (RS) (open circles). The solid lines 
represent the least square fits to the data points. 
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Table 3. Spectral parameters of CC Eri (2T models ). 
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Table 4. Spectral parameters 
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CC Eri (CEM RS-model). 
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Even for the bins excluding the flare, the spectral re- 
sults derived with either the 2T or CEM RS model all show 
the presence of high temperature (~ 10 r K) plasma around 
CC Eri. This is consistent with the previous detections of a 
number of flare stars in quiescence with the EXOSAT ME 
and the Einstein IPC (e.g. Pallavicini et al. 1990; Schmitt 
et al. 1990). 



4 DISCUSSION 

4.1 The 1990 July 10 Event the Counterpart 
of a Solar Two Ribbon Flare? 

Pallavicini et al. (1990) classify the 32 flares observed on 
dMe stars with the EXOSAT LE detector as impulsive flares 
and long decay flares according to their X-ray decay times. 
Those in the first group, which are similar to solar compact 
flares, have rapid rise and decay times. Those in the second 
group, which are similar to solar two ribbon flares, have de- 
cay times of roughly one hour or longer. However the X-ray 
luminosity and the total energy released during the flares 
of both types are larger than their solar analogs by sev- 
eral orders of magnitude. The very energetic stellar flares 
have volumes and inferred magnetic field strengths which 
are much larger than their solar counterparts. 

The R OS A T observations presented in the previous sec- 
tion indicate that the X-ray spectrum of CC Eri is vari- 
able and that the temperature and emission measure of 
the corona plasma to some extent rise and decay with the 
source intensity. Such a correlation could imply that a flare 
is present and the corona plasma is being heated during the 



flare rise and cooling during the decay phase. The event de- 
tected on 1990 July 10 from CC Eri (see Fig. 1), with a one 
hour rise time and a two hour decay time, may be a counter- 
part of a solar two ribbon flare. The energy released during 
the flare in the X-ray passband alone is ~ 2.9 x 10 33 erg, 
after subtracting the post-flare quiescent flux, which is an 
order of magnitude larger than the total energy in a typical 
large solar flare (~ 10 32 erg over ~ 10 4 s, e.g. Priest 1981). 

In solar two ribbon flares an entire arcade of magnetic 
loops is disrupted by an eruptive event created by the com- 
plex thermal and non-thermal processes which occur dur- 
ing the flare rise phase. This open field configuration subse- 
quently closes back, leading to the formation of a growing 
system of magnetic loops whose footpoints are anchored in 
the bright H a ribbons. It is believed that the growing system 
of loops is a consequence of magnetic reconnection (Forbes, 
Malherbe & Priest 1989). The reconnection theory assumes 
that the open field structure relaxes back to a closed con- 
figuration of lower energy and the excess magnetic energy 
released in this merging process appears as the thermal en- 
ergy of the bright X-ray loops (see Forbes et al. 1989 for re- 
cent development in the theory). For example, based on this 
hypothesis an analytical model for the reconnection process 
was developed by Kopp & Poletto (1984) and was applied 
with some success to the decay phase of several solar two 
ribbon flares. Such a time dependent model has also been 
applied to long duration stellar flares observed from the flare 
stars EQ Peg and Prox Cen with the Einstein and EXOSAT 
observatories (Poletto, Pallavicini & Kopp 1988). It was 
shown that the model is capable of reproducing the energy 
release rate and temporal evolution during the flare decay 
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phase. Given a number of assumptions a set of well-defined 
physical parameters for the emitting regions of EQ Peg and 
Prox Cen was derived by fitting the analytical model to the 
data. 

In order to investigate the physical nature of the flare 
on CC Eri, we have adopted the magnetic reconnection the- 
ory of Kopp & Poletto (1984) to model the decay phase of 
the event on 1990 July 10. Other flare models could be cho- 
sen, and here we show only that this particular model leads 
to solutions with physically plausible parameters. We have 
made the same sets of assumptions as in Poletto et al. (1988) 
regarding the location and size of the active region. Details 
of the model and the assumptions made can be found in 
Kopp & Poletto (1984) and in Poletto et al. (1988). Briefly, 
the energy release rate per radian of longitude during the 
magnetic reconnection is given by: 

^ = (l/8«-)2n(n + l)(2n + if R%B 2 m 

„2nr (2n+l) _ ^ 

[^.W/^^. ^^^^^^^ d^) (3) 

where I\ t 2 = f -P^(#)d(cos 8), and P n (9) is the Legendre 
polynomial of degree n, which is 5, 9, 17 and 35, correspond- 
ing to the latitude width of the region 33°, 20°, 10°, and 5° 
respectively. The longitude width is assumed to be 1.5 times 
the latitude width, as typically seen in solar two ribbon 
flares. B m is the maximum surface field within the region 
and is the stellar radius of CC Eri, which is estimated 
to be ~ 0.60i?Q from the absolute visual magnitude given 
by Evans (1959) and Pettersen (1991). The time dependent 
function y is of the form y(t) = I + (i? m /i?*)[l — exp(— t/to)]. 
The parameter H m is the maximum height reached by the 
reconnection point during its upward movement and is as- 
sumed to be equal to the separation of the loop footpoints. 
to is the time constant of the energy release, which is found 
from the modelling and has a mean value of 17 ks. Follow- 
ing Poletto et al. (1988) it is assumed that 10 percent of 
the energy released by the magnetic reconnection process is 
radiated in the X-ray band. 

To estimate the energy release rate during the 1990 July 
10 event, we have fitted the CEM RS model to the difference 
spectra of CC Eri, with Nu fixed at 2.6 x 1 18 cm -2 and a 
at 1 .15. The X-ray light curve shown in Fig. I has been sub- 
divided into 17 time intervals (labelled spectral group 3-19 
in Table 2) over which the PHA spectra have been accumu- 
lated. The difference spectra are obtained by subtracting the 
quiescent spectrum measured in group 19 from the spectra 
3-18. As given in Table 5, in the decay phase of the 1990 
July 10 event, the temperature decreases from 2.8 x I0 r K 
to 1.2 x 10 K and the emission measure from 4.2 x I0 52 to 
0.1 x I0 52 cm" 3 . 

We have fitted equation (3), after integration along the 
longitude and multiplied by a constant 0.1 , to the energy re- 
lease rate measured with the PSPC during the decay phase. 
In Fig. 4 we plot the observed energy release rate together 
with the best fit curve which has little dependence on the 
polynomial degree n. It is obvious that this magnetic re- 
connection model could account for the decay phase of the 
flare-like event on 1990 July 10 except in the period 19:39- 
19:57 where a second bump exists. (This may be caused 
by a second flare). In Table 6, we list a set of geometri- 



Table 5. Parameters of difference spectra (CEM RS-model). 
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cal and physical parameters of the X-ray emitting region 
derived from the model fitting. The maximum surface mag- 
netic field within the region B m ranges from 252 Gauss to 
1474 Gauss, depending on the size of the flare region. This 
is the minimum value required to account for the flare en- 
ergy. T4ise , the initial upward velocity of the reconnection 
point, defined as V vlse = H m /to, is similar to those derived 
for EQ Peg and Prox Cen (Poletto et al. 1988). The electron 
density N e , which is evaluated for the decay phase from the 
observed emission measure (see Table 5) by approximating 
the X-ray emitting loop arcade with a semi-cylinder whose 
volume at time t is ~ 1,5l ^ Hm [I — exp(— t/to)] 2 , decreases by 
about an order of magnitude from the flare peak to the qui- 
escent state. The values of N e and B m required lie between 
those found for the flares on EQ Peg and Prox Cen (Polleto 
et al. 1988), but the values of H m are larger because the 
radius of CC Eri is larger. The last column in Table 6 rep- 
resents a limiting case since it gives a radiation loss (0.2-2 
keV band) that just exceeds the stellar surface flux. 

There are no measurements of the surface magnetic field 
of CC Eri. Measurements for other stars suggest that the 
fields are close to the equipartition values (Saar, Linsky & 
Beckers 1986). Using T e jf as given in Table I, the equipar- 
tition field would be 2630 Gauss. The observations of 17 
stars with measured fields and area filling factors (see Mon- 
tesinos & Jordan 1993) give a relation between log / s and the 
Rossby number Ro (Ro = P ro t/r c , where r c is the turnover 
time at the base of the convection zone, calculated accord- 
ing to Noyes et al. 1984). With Ro = 6.0 x I0 -2 , the filling 
factor is predicted to be / s = 0.7. Thus the minimum values 
of B m in Table 6 lie below the equipartition field, as they do 
for the flare on Prox Cen (Poletto et al. 1988). 

The parameters in Table 6 can be compared with 
those of large solar two-ribbon flares (See Bray et al. 1991; 
Pallavicini et al. 1990). Even the largest solar flares have 
areas which are less than about 0.2% of that of the solar 
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Table 6. Flare parameters of CC Eri. 
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Figure 4. Evolution of the energy release rate of the two ribbon flare. The crosses represent the energy release rate measured by the 
PSPC and the solid line is the best fit curve of the magnetic reconnection model. 



hemisphere. With an area of less than 1% the main differ- 
ence in the flare on CC Eri would be the high electron pres- 
sure implied (greater than 2 x 10 19 cm -3 K, compared with 
values an order of magnitude lower in a solar flare). Alterna- 
tively, a pressure similar to that in a solar flare is found only 
with a substantially larger area (8% of the hemisphere) and 
larger height, than typical of a large solar flare. Thus in this 
model, which has continual energy input, the flare on CC Eri 
has either a higher density or a larger volume than a large 
solar two-ribbon flare. From this model alone any solution 
with an area between 1 and 10% of the stellar hemisphere 
gives physically possible parameters. In Section 5.1 we dis- 
cuss an alternative model with no heating beyond the flare 
maximum, in attempt to restrict the possible solutions. 



4.2 The 1990 July 10 Event Modulation of 
Stellar Rotation? 

As for other BY Dra variables, the optical photometric and 
spectroscopic peculiarities of CC Eri are attributed to spots 
on the star surface (Krzeminski 1969; Bopp & Evans 1973; 
Busko et al. 1977). The cause of the variation is the presence 



on the surface of magnetic starspots (in active regions) with 
a temperature a few hundred degrees cooler than the sur- 
rounding photosphere. It is expected that the starspots give 
rise to the optical light minimum and the associated active 
regions to the X-ray flux maximum. The starspot model, 
via a solar analogy, may provide a ready explanation for 
the X-ray activity observed with the PSPC. The X-ray flux 
increases and decreases when the active regions associated 
with the spots pass in and out of the line of sight. Based on 
this scenario we have created a starspot model to reproduce 
the X-ray light curve as observed with PSPC on 1990 July 
10. Since the spot areas are relatively large, we assume that 
the active regions have the same size as the spots. In prac- 
tice, this may underestimate the size of the active regions. 

In a spherical coordinate system, 7, the angle between 
the line of sight and the normal line of an X-ray emitting 
area cell can be written as 

cos 7 = cos(ip — (f>) sin 8 sin i + cos 8 cos i (4) 

where <f> is the stellar longitude and 8 the co-latitude; 
if) is defined as the rotational phase and i as the orbital 
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0.08 0.1 
Stellar Rotation Phase 

Figure 5. X-ray light curve of CC Eri as a function of the stellar rotation phase. The synthetic light curve from the starspot model is 
displayed as a solid line. The contributions of the north and south spots are shown by the dash-dot-dash line and by the long dash line 
respectively. 
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Table 7. Spot parameters of CC Eri. 
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to —51.5°. Note that the fluxes required for the central part of the south spot exceed the stellar surface flux; see text 
for discussion. 



inclination which is 42° (Evans 1959; Bopp & Evans 1973). 
The X-ray region is visible if the angle 7 is less than 90°. 

Two starspots (which may be groups of spots), one lo- 
cated in the north hemisphere and another in the south, are 
required to reproduce the X-ray light curve on 1990 July 
10. The north hemisphere spot has a uniform flux distribu- 
tion, and the south spot has a flux distribution in which 
the flux is high in the central region (<f> ~ 18° — 28°, and 
8 ~ 131.5° — 141.5°), low in the surrounding area. This re- 
gion of concentrated flux is required to match the observed 
X-ray variation. In Table 7 we list the parameters of these 
spots. Each spot is divided into 100x100 X-ray emitting cells 
and the X-ray light curve is synthesized by adding the con- 
tribution of each cell according to equation (4). The shape 



of the X-ray light curve is largely determined by the location 
and size of the spot in latitude. The position of the X-ray 
flux maximum and minimum in the stellar rotation phase is 
a function of the longitude. We have adopted an ephemeris 
from Evans (1959) in which the stellar rotation period is 
1.56145 days and the phase is at JD2430000.0. The peak of 
the observed X-ray flux occurs at rotation phase ip = 0.06. 
The longitude location of the spots is thus determined by 
maximizing equation (4). The size of the northern spot and 
its latitude location are not independently known and the 
three values used (see Table 7) are simply adopted from 
Bopp & Evans (1973) who modelled the optical light curves, 
over three periods of time, with these parameters. We have 
generated optical light curves using the same limb-darkening 
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coefficients and ratio of spot temperature to photospheric 
temperatures as those used by Bopp & Evans (1973). The 
predicted ratios (see Table 7) of the optical minimum to 
maximum luminosities L mtn / 1 Lmax are consistent with the 
optical observations. 

The X-ray light curve of CC Eri shown in Fig. 1 has 
been folded with the stellar rotation phase and is plotted in 
Fig. 5 together with the synthetic light curve and its two 
components. In the rotation phase 0-0.17, the size and loca- 
tion of the north spot has no effect on the shape of the X-ray 
light curve but gives a constant flux. The X-ray variation in 
this phase interval is completely determined by the loca- 
tion, size and the X-ray flux of the south spot. By contrast, 
the optical photometric variation is determined largely by 
the north spot. Although it is possible to choose parameters 
for the south spot that can reproduce the X-ray variations, 
as shown in Fig. 5, the solutions for the brightest region 
are not physically acceptable since even the radiation losses 
(0.2-2 keV band) alone exceed the stellar surface flux. Thus 
it seems unlikely that the variations observed are due to ro- 
tational modulation of activity related to starspot regions. 



5 MODELLING FROM THE EMISSION 
MEASURES AND TEMPERATURES 

Since emission measures and temperatures have been de- 
rived, some further modelling is possible. First, the change 
in the emission measure and temperature during the decay 
phase of the flare can be used to deduce geometrical param- 
eters and limits to energy fluxes, under the assumption of 
no heating. This is an alternative model to that discussed in 
Section 4.1, which assumed continued energy input through 
magnetic reconnection. Secondly, the spectral groups before 
and after the flare can be used to investigate the average 
coronal conditions. The emission measure distribution de- 
rived from an energy balance model can be used to make 
comparisons with the value found at 1 — 2 x 10 5 K found 
from analyzing spectra obtained with the IUE, which are 
now in the public archive. 



we use only the changes in Em(V) and T e , and the aver- 
age values, between (a) groups 8 to 11 and (b) 13 to 15, to 
examine the possible values of L and V. 

If F c (T e ) ^> i*R,(Te), then L can be found, giving values 
of (a) ~ 3.1 x 10 10 cm and (b) ~ 1.8 x 10 10 cm. For com- 
parison, the pressure-squared isothermal scale height at the 
mean temperatures are ~ 4.0 x 10 10 cm and ~ 2.8 x 10 10 
cm. Only the product N e A 1/2 (A - area of one footpoint) 
can then be found and this has values of (a) 6.4 x 10 20 cm -2 
and (b) 5.5 x 10 20 cm" 2 . 

If -Fi),(Te) ^> F c (T e ), then the volume can be found and 
the values are (a) 2V = 1.5 x 10 30 cm 3 and (b) 1.1 x 10 30 
cm 3 . (These can be compared with the volume of a spher- 
ically symmetric corona, of thickness corresponding to the 
pressure-squared isothermal scale height at the average qui- 
escent temperature given in Section 5.2, i.e. ~ 2 x 10 32 cm 3 ). 
In this case the density can also be found and the values are 
(a) N e ~ 1.3 x 10 11 cm" 3 and (b) 1.1 x 10 11 cm" 3 , giving 
pressures comparable with those in a large solar flare. 

Increasing the radiative term relative to the conductive 
term leads to larger values of L and V, but lower values of 
NeA 1 / 2 and N e . Adding a heating term would, for the same 
values of dT/dt, lead to smaller values of L and V. 

Assuming that F c = _Fr the total energy losses are (a) 
5.3 x 10 9 erg cm" 2 s" 1 and (b) 2.6 x 10 9 erg cm" 2 s" 1 , about 
a factor of three to five lower than the stellar surface flux 
(1.3 x 10 10 erg cm" 2 s" 1 ). 

Although the assumption of no heating may not be ap- 
propriate early in the decay, it is of interest to compare the 
parameters for groups 13 to 15, late in the flare, with those 
in the heated model, since in the latter case the heating 
decreases exponentially according to equation (3). Table 8 
gives the parameters derived in the two models, for the case 
of F c > _Fr since in the heated model it is assumed that only 
10% of the total energy loss is in the form of radiation in 
the X-ray band. From Table 8 it can be seen that the two 
models agree best for the reconnection model which has an 
area of 8%. This area is within the range of star spot areas 
deduced by Bopp & Evans (1973), as given in Table 7. 



5.1 The Time- Varying Component 

If the variation in the flux and temperature are interpreted 
as a flare, the difference spectra can be used to estimate ei- 
ther the length or the volume of the region. The full equation 
for the cooling rate, including a varying emission measure 
and temperature has been given by Veck et al. (1984). The 
conductive flux can be estimated by balancing it against 
the radiation below the region of maximum temperature. 
The equation for the change in T e then becomes 

dT e _ -7.2 x lO 4 ^ 1 25 /L - (0A9Em(V) 1/2 /V 1/2 T e ) 
dt ~ 1 - (dlog(Em(V)/V)/3 dlogTe) ^' 

where Em(V) is the volume emission measure for half 
the loop, L is the loop half-length and V is the half-loop vol- 
ume. The denominator on the R.H.S. of equation (5) must be 
> 0. Using the data given in Table 5, dlog Em(V)/3 dlog T e 
satisfies this condition, between groups 7 to 11, and between 
groups 13 to 15. However, any increase in V must be small, 
otherwise the denominator becomes < 0. Since the uncer- 
tainties in the individual values of Em(V) and T e are large, 



5.2 Average Conditions 

The mean coronal emission measure and temperature are 
derived from the ten spectral groups, 1 to 4 and 15 to 20. 
(Using a more restricted set 1 to 3 and 17 to 20 makes no 
significant difference). Studies of other stars with ROSAT 
and EUVE (Jordan et al. in preparation) suggest that of the 
various spectral fits, the higher temperature of a two tem- 
perature fit has the most physical significance. We adopt the 
results from the RS code. The mean values of the apparent 
emission measure and temperature are: 

Em(V) = J N 2 dV = 1.06 x 10 52 cm" 3 (6) 

and 

T e = 9.46 x 10 6 K (7) 

In practice, only some fraction of the emitting volume 
is observed, since the star shadows the rest. The fraction of 
a spherically symmetric corona that is observable is given 
by 
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Table 8. Comparison of parameters derived from the model with no heating and the reconnection model. 





Model with 


no heating 


Area: 


Reconnection model 

8% 


3% 




L 


N e A 1/2 


L 


N A 1 / 2 L 


N e A 1/2 




(10 10 cm) 


(10 20 cm" 2 ) 


(10 10 cm) 


(lO 2 " cm" 2 ) (10 10 cm) 


(10 20 cm" 


Groups 8-11 


3.1 


6.4 


1.4 


1.3 0.84 


1.7 


Groups 13-15 


1.8 


5.5 


2.0 


7.4 1.2 


9.5 



G(r) = 0.5(1 + (l-(^/r) 2 ) 1/2 ) 



(8) 



Allowing also for the radial extent of the corona, the 
apparent emission measure as a function of r is 



Em(r] 



App 



Em(r 



The "true" emission measure is then 
Em(V) 



^R 2 G(r)f(r) 



(9) 



(10) 



where 



/W = (^) 2 (ii) 

In a plane parallel approximation to a uniform, non- 
extended corona, the value of G(r) is 0.5, and f(r) is 1.0, so 
that the maximum value of Em(r)T is 9.67 x 10 29 cm -5 . 

The electron pressure can then be found by writing 
Em(r)T as 

Em(r) T = P 2 |-T- 2 (12) 

where P e is in cm -3 K and H /2 is the pressure-squared 
isothermal scale height, given by 

y = 7.09 x lO 7 ^* 1 (13) 

This simple approximation gives the maximum pressure 
and density in a plane parallel corona, since including the ra- 
dial terms leads to lower values. (The pressures and densities 
would of course be even higher if the emitting area was re- 
stricted to only part of the corona, but there is not enough 
information to constrain such models). 

In the plane parallel approximation the observed values 
of Em(V) and T e lead to P e = 6.8 x 10 16 cm" 3 K and N e = 
7.2 x 10 9 cm -3 . However, this approximation will not be 
entirely appropriate because the isothermal scale height is 
significant compared with the stellar radius. The average 
radius of the emitting corona can be estimated from H/2. 
(The height of the base of the corona is likely to be small, 
by analogy with the Sun). With H/2 = 1.89 x 10 10 cm, one 
finds (r/R-t) < 1.45 and G(r) < 0.86. Thus the lower limits 
to the true emission measure and the electron pressure are 
Em(r) T > 2.7 x 10 29 cm -5 , and P e > 3.6 x 10 16 cm -3 K. 

The emission measure distribution below the coronal 
temperature is not known. It is possible to calculate the dis- 
tribution by assuming that there is no non-thermal energy 
deposited below the corona, so that the radiation losses are 
balanced by the net conductive flux. In the spherically sym- 
metric case the energy balance gives 

r 2 ar Li 
where 



-0.87^fP rad (T e 



(14) 



Fc(T e . 



-K T e 



5/2 



dT e 
dr 



(15) 



and K = 1.1 x 10 -6 erg cm -1 K _r/2 s _1 . The form of 
radiative power loss coefficient adopted is T ra d(T e ) = e/T e = 
2.33 x 10 _16 /Te erg cm s , to reproduce the RS code at 
10 r K. Substituting 



dTe/dr 



P 2 



T e yiEm(0.3) 



(16) 



where, Tm(0.3) is the average Em(r)T over AlogT e = 
0.3 dex, to match the typical region of line formation, leads 
to 

V2£m(0-3) ] ^ T 2 (17) 



d[(^) 2 «oT 3 ' 2 



Expressing Tm(0.3) in terms of the apparent emission 
measure, Em(r), equation (17) can be written as 



dlog(Tm(r)/G(r)) 
d log T e 



d log r 3 d log P e 
4 1 h 2 s—l 

d log T e 2 d log T e 

_2 £q ^ 4 Err^ 1 
ko r T e 5/2 G(r) T e 2 

Our computer programme uses the apparent coronal 
emission measure, temperature and pressure to calculate 
Em(r) in an iterative manner, making use of the subsidiary 
equation (16) and the equation of hydrostatic equilibrium, 
expressed as 

d_P e fim p g^\/2 Em(r) 1 
d% = 



(19) 



k G(r) f(rf 

The radial factors are removed in the plane parallel ver- 
sion, and G(r) is set to 1/2. The absolute height scale is 
initially arbitrary. 

The resulting emission measures are then used together 
with equations (16) and (19), to calculate the run of pres- 
sure and temperature with height above the chromosphere, 
using also, below T e = 2 x 10 5 K, emission measures de- 
rived from LUE spectra. Here we are concerned only with 
a comparison of the observed and calculated values of the 
emission measure at T e = 1 — 2 x 10 5 K. The full emission 
measure distribution and models of the chromosphere to the 
corona will be discussed in a later paper. The relevant emis- 
sion measures are those found from the lines of C IV and 
N V. We have extracted 12 low resolution short wavelength 
IUE spectra from the archives to find the average surface 
fluxes. (See also Byrne et al. 1992). The emission measures 
were found as in Jordan et al. (1987), using the "solar" set 
of ion abundances from Jordan (1969), element abundances 
from Grevesse, Noels & Sauval (1992), the same collision 
strength for N V, but that for C IV reduced by a factor of 
1.17. The resulting values of Em(0.3) are 3.2 x 10 28 cm" 5 
for C IV and < 6.6 x 10 28 cm -5 for N V. 
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The plane parallel models are very sensitive to the coro- 
nal values of Em(r) and e adopted when Em(r) is near a 
critical maximum value. Adopting Em(r) = 9.7 xlO 29 cm -3 , 
does not give an energy balance solution, but a slight reduc- 
tion to 8.2 x 10 29 cm -5 gives a solution with an emission 
measure of 3.6 x 10 28 cm -5 at 2 x 10 5 K, similar to that 
derived from the IUE spectra. A small reduction in e would 
achieve the same result. The coronal pressure is 6.4 x 10 16 
cm -3 K and the coronal temperature (rounded to 10 r ) is 
reached by a height of (h + R*) = 5.1 xlO 10 cm (r / R* = 1.2). 

In the spherically symmetric model the starting appar- 
ent emission measure is defined by equation (9), and the geo- 
metrical factors are then computed as part of the iterations. 
Because of the factor f(r) (= 1.47) the true emission mea- 
sure is smaller (2.6 x 10 29 cm -5 ), leading to a lower coronal 
electron pressure (3.8 x 10 16 cm -3 K), but the temperature 
gradients (from equation (16)) are similar. The model then 
has the same radial extent, but the emission measure pro- 
duced at T e ~ 2 x 10 5 K is two orders of magnitude lower 
than that observed. 

The plane parallel model is clearly not entirely self- 
consistent since the radial extent should be taken into ac- 
count. But neither is spherically symmetric model satisfac- 
tory. This suggests that more complex geometries should be 
explored, but would require assumptions regarding area fac- 
tors and the magnetic field geometry. However, the present 
models do limit the range of the coronal pressure to within 
about a factor of two. The values found are comparable with 
those in a well developed solar active region and are about 
two orders of magnitude lower than found in the CC Eri flare 
models. Observations with EUVE could in principal further 
constrain the emission measure distribution and thus the 
geometry. 

The observed parameters can be compared with the pre- 
dictions of Hearn's (1975) minimum energy loss (m.e.l.) hy- 
pothesis. The equations used are given in Montesinos & Jor- 
dan (1993) (equations 4.8 and 4.9). With T e = 9.46 x 10 6 K, 
the predicted emission measure and pressure are Em(r) = 
4.8 x 10 29 cm" 5 and P e = 4.7 x 10 16 cm" 3 K. These val- 
ues lie between the 'observed' values with the plane paral- 
lel and spherically symmetric interpretations. As found by 
Montesinos & Jordan (1993) for a larger sample of main- 
sequence stars, the m.e.l. hypothesis gives a remarkably good 
prediction of average coronal properties. 

The observed emission measure and temperature can 
also be used to estimate the energy lost from the corona by 
radiation and conduction. Using equation (4.18) and (4.21, 
or 4.22) from Montesinos & Jordan (1993), these are, respec- 
tively, 2.1 x 10 r erg cm" 2 s" 1 and 1.0 x 10 8 erg cm" 2 s" 1 , 
in the plane parallel interpretation. The total energy loss 
follows the trend with the Rossby number and g*, found by 
Montesinos & Jordan (1993) (equations 4.30). The radiation 
flux would be smaller in the spherically symmetric approxi- 
mation. Thus even in the "average" corona, the energy loss 
are ~ 1% of the stellar surface flux, significantly larger than 
in the Sun. 

Using the definition of the plasma /?, the coronal mag- 
netic field can be found from the emission measure and tem- 
perature (see equation (4.4) in Montesinos & Jordan 1993), 
giving B c = 21/?"° 5 Gauss (plane parallel model). 



6 CONCLUSIONS 

CC Eri was observed in the periods 1990 July 9-11 and 1992 
January 26-27 with the PSPC detector on board the ROSAT 
satellite. These high quality data give the first information 
on the temporal and spectral variability of CC Eri in the 
X-ray energy band. 

During the R OS AT observations the X-ray intensity of 
the source is variable on timescales from a few minutes to 
hours. The X-ray luminosity is in the range 2.5 — 6.8 x 10 29 
erg s" 1 , which is similar to values found from previous mea- 
surements with the EinsteinlPC and EXOSATLE. On 1990 
July 10 an X-ray flare-like event was detected with an ex- 
ponential rise time of about one hour and a decay time of 
about two hours. 

The X-ray spectrum of CC Eri can be well reproduced 
with either 2T or CEM (RS) models and the spectral results 
derived with either show the presence of high temperature 
(~ 10 r K) plasma around CC Eri. We find that the X-ray 
spectrum is also variable and the variations of the emission 
measure and to a lesser extent, the temperature, are corre- 
lated with the source intensity. 

The X-ray variability of CC Eri observed with the PSPC 
may in principle be caused by flaring events and/or by the 
rotational modulation of active regions. However, the pa- 
rameters required for the region causing X-ray variations are 
physically unrealistic. Thus it seems more likely that some 
sort of flare was observed. A model developed for solar two 
ribbon flares in which heating is produced by magnetic re- 
connection requires either that the flare has a larger volume 
and percentage area than a solar two ribbon flare, or higher 
densities. A comparison between the reconnection model pa- 
rameters and an unheated model, towards the end of the 
flare, suggests that the solution with a larger area (~ 8%) 
and volume, and electron pressures similar to those in solar 
flares is the more likely. The larger flare area on CC Eri, 
compared with the Sun, may not be surprising given the 
larger optical star spot area. Thus a two ribbon flare pro- 
vides one possible explanation of the variations observed. 
The occurrence of a flare would, by analogy with the Sun, 
require the presence of an active region, but the contribu- 
tion of such a region to the 'quiescent' emission cannot be 
determined. If the non-flaring emission is attributed to an 
'average' corona, then the electron pressure and emission 
measure are similar to those predicted by the minimum en- 
ergy loss hypothesis, which gives an adequate description 
of a range of main-sequence stellar coronae. The electron 
pressure in the average corona is ~ 4 — 6 x 10 16 cm" 3 K, 
comparable with that in a well-developed solar active region. 
The energy flux (radiation plus conduction) is ~ 1% of the 
stellar surface flux, substantially larger than for the Sun. 

Ideally, CC Eri needs to be monitored simultaneously 
in optical, ultra-violet, and X-ray wavelengths over a com- 
plete binary cycle in future observations. Such observations 
could be used to study both flares and rotational modula- 
tion. Multi-wavelength observations could be used to map 
the X-ray distribution over the stellar disc and constrain 
the area covered by active regions. Observations with EUVE 
could be used to derive information on the emission measure 
distribution and hence place limits on the geometry of the 
emitting regions. 
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